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ABSTRACT 



C/^ . Received: August 2009 / Accepted 

^' : 

O Context. Magellanic Clouds are of extreme importance to the study of the star formation process in low metallicity environments. 

, Aims. In this paper we report on the discovery of pre-main sequence candidates and young embedded stellar objects in Nil 
^ ' located in the Large Magellanic Cloud to cast light on the star formation scenario. We would like to remind that this comparison 
I I ' is complicated by the presence of a large age dispersion detected in the fields. 

. Methods. Deep archive HST/ACS photometry is used to derive color-magnitude diagrams of the associations in N 11 and of 
' the foreground field population. These data are complemented by archive IR Spitzer data which allow the detection of young 
^ [ embedded stellar objects. The spatial distribution of the pre-main sequence candidates and young embedded stellar objects is 
compared with literature data observed at different wavelengths, such as Hq and CO maps, and with the distribution of OB and 
Herbig Ae/Be stars. The degree of clustering is derived using the Minimal Spanning Tree method and the two point correlation 
function to get insights about the formation process. 
. Results. A large population of pre-main sequence candidates is found in N 11. Their masses are in the range of 1.3—2 Mq 
for ages from 2 to 10 Myr. Young embedded stellar objects having ages of 0.1 — 1 Myr are found to be intermixed with the 
candidate pre-main sequence stars. The spatial distribution of the stars shows that this region is the product of clustered star 
formation. No significant difference is found in the clustering degree of young blue main sequence stars and faint pre-main 
sequence candidates, suggesting that they might be part of the same formation process. 

Conclusions. The data suggest that the star formation in the region is a long-lasting process where stars from 0.1 to 10 Myr are 
widely distributed. 
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1. Introduction 

Several plausible scenarios are presented in the literature 



to explain how i solate d stars form (e.g lBonnell et al.l . l2004 : 
2OO3I) . A detailed understanding of how 
converted into stars is still miss 



McKee fc Tan 
a molecular cloud is 



ing (see among others lElmegreen . l2000t iDib erail l2007l : 



Krumholz et al.l . l2007f ). Interstellar turbulence is found to 



be very efficient in sweeping up molecular gas-forming 
massive structures which in turn can undergo a large-scale 
gravitational collapse. Supersonic turbulence can counter- 
balance gravity on a global scale, but provokes c ollapse 
locally (|Sasaol Il973t iKlessenl . l2003l : iLarsonl . |2003|) . As a 
consequence, clusters are formed in a hierarchical fash- 
ion with subclusters which eventually merge to build up 
the final condensed object ( Schmeia et al. . 2009f) . Stellar 
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feedback is expected to strongly influence the process. 
Young massive stars inject energy into the nearby inter- 
stellar medium, heating and compressing the surrounding 
gas. This process can have destructive or constructive ef- 
fects, depending on the balance between heating and grav- 
ity, but it is still not clear what regulates it. The pres- 
ence or absence of turbulent feedback directly relates to 
the physical mechanism of star formation and determines 
whether stars are generated by t he formation and col- 



lapse of discrete prot o stella r cores (jKrumholz et al.l . 12005 



Padoan fc Nordlund . 
(|Bonnell et allbooj . 



20021 ) or by competitive accretion 



To analyze entire star-forming regions including dif- 
ferent components (molecular clouds, ionized gas, stars.. ) 
can cast light on the formation scenario. Clustering prop- 
erties of the young population can give information about 
the structure of the interstellar medium from which it 
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formed. The Magellanic Clouds are ideal laboratories to 
study the process of star formation in detail, due to their 
proximity. In the Large Magellanic Cloud (LMC) several 
regions of active star formation are found. Nil is the 
second largest nebula of the LMC after the 30 Doradus 
Nebula. It is located at the North- Western corner of the 
LMC. 

This region is often presente d as one of the best exam- 
ples of trigger e d star formation dWa lborn fc Pa rkeii 1992 ; 



Rosado et al 
20071) . 



199d iHatano et all 2006; Mo kiem et al 
The basic idea is that the association LH 9 has 
triggered star formation in LH 10 and LH 13 (see discus- 
sion in the following). 

Here we discuss the stellar content in the region of N 1 1 
using HST ACS/WFC archive data. These observations 
are complemented by Spitzer archive data: while HST ob- 
servations can give information about faint, exposed pre- 
main sequence candidates (PMSs), IR data allow to detect 
embedded young stellar objects(YSOs). Regions of active 
star formation can be detected, pointing out the differ- 
ences between the areas where the nebulae are located 
and the surrounding fields. 

This paper is part of a project aimed to cast light on 
the process of field and cluster star formation in the 



Magell anic Clouds. We quote among others iChiosi et al.l 
(|2006a[) where the cluster and field star formation in the 
central part of the Small Magellanic Cloud (SMC) is stud- 
ied and the correlation between young objects and their in- 
terste llar environment is discussed. In Chiosi fc Vallenaril 
(j2007f ) the star formation rate in the SMC is found to be 
quiescent at ages older than 6 Gyr. 

This paper is organized as follows: in Sect. [2] the clus- 
ter and association content of the area is presented; in 
Sect. El the HST ACS/WFC and Spitzer archive data are 
presented; in Sect |4] the methods used to analyze the re- 
gion are discussed; in Sect [5] the field population is studied 
and the star formation history is derived. In Sect [6] an es- 
timate of the interstellar extinction is given; in Sect [7] the 
CMDs of the most relevant associations and clusters are 
presented and their stellar contents and ages are discussed; 
in Sect[H]the degree of clustering of the stellar populations 
at different ages is derived. Concluding remarks are drawn 
in Sect[9l 

2. Presenting the region: clusters and associations 
in Nil 

This section summarizes literature information about the 
cluster and associations in Nil. The complex has a diam- 
eter of 45', corresponding to a linear ext ent of 655 pc, if 
the d istance modulus (M — m)o = 18.5 (IKeller fc Wood 



200i). This region has a peculiar morphology: a central 
hole with no emission is surrounded by several bright 
nebulae and filaments. These structures are observed at 
the Ha and fOHI] 5007 A wavelength s and in the radio 
(|Rosado et al1 . ll996l : lls"rael et al.L 12003^. 
The different emission areas are classified according to the 
emission brightness starting from Nil A to B,C,D,E,F and 




Fig. 1. The region containing the studied HST/ ACS fields 
(boxes) is schematically represented. The ellipses indicate 
objects of particular interest, as clusters and associations. 
The numbers refer to the object ID listed in Table [T] 

I. The central cavity has been evacuated by the OB stars 
in the association LH 9, located near the center of Nil. 
The OB association LH 10 in NllB is to the north of LH 9, 
while LH 13 in NllC is to the east of LH 9. LH 14 in NllE 
is situated to the northeast of LH 10 and is outside the 
studied region. 

The Bica et al. I ()l999l) Catalog lists 49 clusters and as- 
sociations in the whole region of Nil. This catalog collects 
and cross-identifies stellar and non-stellar objects (nebu- 
lae, supernova remnants) out of more than 30 catalogs. 
In our observed region, 13 objects are detected. Table [T] 
lists the associatio ns and clusters in the observed region 
from the catalog of iBica et al. 1 (Il999l) together with their 
minimum and maximum diameter, and the Bica et al 



( 1999i) classification. For clarity, we recall it. The sequence 
C (Clusters) to A (Associations) refers to the density, 
where high density objects are classified as C, and low 
density objects as A. Objects NA and NC are stellar sys- 
tems clearly related to emission. Figure [1] presents a map 
with the loca t ion of clusters and associations of interest. 
Hatano et al. ( 20061 ) find 127 Herbig Ae/Be star candi- 
dates from near-infrared photometry, mainly in the pe- 
riphery of LH 9. Herbig Ae/Be star are intermediate mass 
prc-main sequence stars (3—7 Mq) having an age range of 
1-3 Myr. 

3. The data 

Here we first present the HST ACS/WFC archive data 
we use together with the data reduction and calibration. 
Then we discuss the infrared Spitzer archive data. 
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Table 1. Clusters and associations of particular interest in the region under investigation. Column 1 gives the iden- 
tifier(ID), Cols. 2 and 3 the coordinates (J2000), Col. 4 gives the classification by Bica et al, Cols. 5 and 6 give the 
minimum and the maximum values of the diameter, Cols. 7 and 8 list the name and the location. 



ID 


Ra 




Deg 






Class 




D 




D 


Name 


Location 


















(arcm 


n) 


(arcm 


n) 














Max 


Min 






1 


04 


56 


38 


-66 


30 


26 


NA 


7 


30 


5 


90 


LH 9 


LMC-Nll 


2 


04 


56 


43 


-66 


31 


28 


NA 


5 


10 


2 


10 


NGC 1760, 

LMC-NllF, 

ES085EN19 


LH 9 


3 


04 


56 


39 


-66 


29 


00 


A 


4 


20 


3 


00 


NGC 1761, 
SL 122, 
ES085SC18 


LH 9 


4 


04 


56 


34 


-66 


28 


25 


C 


1 


00 





90 


HD 32228, 
KMHK 307 


NGC 1761 


5 


04 


56 


23 


-66 


28 


00 


NA 





95 





75 


BSDL 270 


NGC 1761 


6 


04 


56 


30 


-66 


26 


40 


NA 


26 


00 


23 


00 


LMC-Nll, 
LMC-DEM34 


LMC-NIO 


7 


04 


57 


44 


-66 


28 


30 


NC 





30 





25 


HNTl 


NGC 1769 


8 


04 


56 


19 


-66 


27 


19 


A 


1 


30 





90 


BCDSP 1 


LH 9 


9 


04 


57 


45 


-66 


27 


52 


NA 


4 


10 


2 


60 


NGC 1769, LH 13, 

LMC-NllC 

ES085SC23 


LMC-Nll 


10 


04 


56 


09 


-66 


31 


25 


NA 


1 


50 





90 


BSDL 264 


LH9 


11 


04 


57 


17 


-66 


24 


56 


NA 





80 





70 


HT3 


LMC-Nll 


12 


04 


57 


44 


-66 


26 


27 


NA 





65 





55 


BSDL 324 


NGC 1769 


13 


04 


56 


49 


-66 


24 


23 


NA 


5 


20 


3 


60 


NGC 1763, 
IC 2115, 
LMC-NllB, 
SL 125, LH 10 


LMC-Nll 



3.1. HST ACS/WFC photometry of the stellar 
population 

Data for six fields in Nil were taken from the ACS/WFI 
HST archive (PI Mai'z-Apellaniz, Proposal ID=9419). The 
studied fields are partially covering the associations LH 9, 
LH 10, and LH 13. In Fig. [1] we sketch the position of the 
ACS fields in the Nil complex. The control field to sample 
the background and foreground population is located out- 
side the Nil boundaries and is not plotted. Observations 
are made in the filters F435W and F814W with long expo- 
sures of 496 s and short exposures of 1—2 s. We make use 
of drizzled files which already account for image distor- 
tion. The data are then red uced by means of the packages 
DAOPHOT/ALLSTAR bv lStetsonI (|l994h . 

The c orrection for charge t ransfer efficiency is applied 



following ISirianni et al.l (j2005l ). The instrumental magni- 
tudes are transformed to the ACS-V ega system using the 
zero points derived bv lBohhnl (l2007f l. 

Figure [3] presents the CMD of the background popu- 
lation. 

Photometric errors. Photometric errors are es- 
timated in the usual way by means of artificial star 
experiments, which consists of the injection of a large 
number of artificial stars of known magnitude into the 
original image: stars are then recovered through the 
whole reduction pipeline, and recovered magnitudes are 
compared to the original ones. The mean magnitude 
difference between injected and recovered stars, which is 
taken as an estimate of the mean photometric error, is 
plotted in Fig. [2] as a function of the magnitude. 

Completeness factor. The completeness factors A 
are calculated for each field from crowding tests as the 
ratio between the recovered stars in a given magnitude 



F435W 
F814W 



Fig. 2. Photometric errors in the two pass-bands as a 
function of the magnitude, as recovered from artificial star 
experiments. 



interval and the number of original stars in the same in- 
terval. 

In the field population region the completeness is 90% 
at F435W~ 26.0 and F814W~ 25.0. In the associations, 
the mean completeness is better than 90% for F814W< 
24.5 and F435W < 25.0 mag, it becomes 80% at F814W~ 
25.0 and F435W ~ 25.5 mag, and then steeply declines. 
At a changing field, the completeness varies by about 7% 
around the mean value, being slightly worse for the field 
including LH 13. In the following, the analysis of the lumi- 
nosity functions (LFs) is limited to F435W < 25.5 mag (or 
F435Wo < 24.5 mag), where the completeness correction 
is quite negligible. 
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Fig. 3. CMD of the field population where isochrones of 
different ages and chemical composition are superposed, 
namely Age=0.1 Gyr, Z=0.019; Age=0.6, Gyr Z=0.008; 
Age=2 Gyr, Z=0.004; Age=10.0 Z=0.001 Gyr(upper 
panel). In the lower panel, the SFR in the field popu- 
lation is presented as a function of age (in yr) (see text 
for details). 



3.2. Infrared Spitzer data 

The region was observed by the infrared survey SAGE 
(Surveying the Agents of a Galaxy's Evolution) using the 
Infrared Array Camera IRAC (3.6, 4.5, 5.8 and 8 /j,m) and 
MIPS (24 iim) instruments on board of the Spitzer Space 
Telescope (jMeixner et all . l2006l lioOTl ) . AU relevant infor- 
mation concerning those arch ive data, data reduction, an d 
the catalog can be found in iMeixner et al.l ( 20061 2007 ). 
Infrared CMDs and color-color plots will be discussed in 
the following sections. 



4. The methods 

This section discusses the methods used to study the stel- 
lar population in N 11. First, wc summarize the Downhill 
Simplex procedure used to derive the star formation his- 
tory (SFR) of the field population. Second, we describe 
how we subtract the field star contamination from HST 
CMDs of the associations and how we derive the extinc- 
tion. Finally, we discuss how the two-point correlation 
function is used to derive the degree of clustering of the 
stellar population. 

4.1. Field population SFR: Downhill Simplex method 

A detailed description of the use of the Downhill Simplex 
method to derive the SFH of the field population can b e 
found in lChiosi etall (|2006bl) . [Chiosi &: Vallemril (|2007l ). 
Here we recall a few points: 

(1) To infer the SFR, theoretical CMDs in different 
age ranges are simulated. The simulations include the 
spread which is due to the observational photometric er- 
rors, reddening, and the effect of photometric incomplete- 
ness. In each age bin, 10000 to 15000 synthetic stars are 
generated down to the photometric completeness limit 
set at the magnitude where the incompleteness factor 
is 50%. The synthetic CMDs stand on the sets of stel- 
lar tracks, is o chron es and single stellar populations by 
Girardi et al. I (|2002h . We assume the chemical enrichment 
history in agreemen t with 'Pa gel fc Tautvaisiene ( 19991 ). 
Carrera et al.l (|2008l ). Table [2] lists the ages and metallici- 
ties of the synthetic stellar populations used in the simu- 
lations. 

(2) The extinction along the line of sight of the field 
population is estimated by means of the fit of the CMD 
with a set of synthetic stellar populations. 

(3) The relative contribution of different popula- 
tions to the total CMD and luminosity/color functions 
to be compared with the observational ones, in other 
words the SFR, is derived from a minimization algo- 
rithm which makes use o f the Downhill Simplex method 
(jHarris fc ZaritskvL liool ) . 

(4) To provide the required constraints to the mini- 
mization procedure, we split the observational CMD into 
a number of suitable magnitude-color bins. Recent work 
concerning the determination of the SFR from the CMDs 
has pointed out the importance of using a binning that 
takes the various stellar evolutionary phases into ac- 
count, as w e U as t he uncertainties of t he stellar models 



(|Rizzi et all . 120021 : IChiosi et all l2006al ). To avoid spuri 



ous solutions, the fit is mainly based on the main sequence 
magnitude. To properly sample the age of the population, 
the magnitude bin is variable, being 0.5 mag from 16 to 
19 mag, and 0.25 mag from 19 mag to 25 mag. 

4.2. Field star subtraction 

Field subtraction is a critical issue to isolate the popula- 
tion of N 11. Field stars are statistically subtracted from 
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Table 2. Ages and metallicities of the synthetic popula- 
tions in use. 



Age 


(Gyr) 




Z 






0.08 


- 0.12 


0.006 







010 


0.12 


" 0.30 


0.006 


— 





010 


0.30 


- 0.40 


0.006 


— 





010 


0.40 


— 0.50 


0.006 







010 


0.50 


- 0.60 


0.006 







010 


0.60 


- 0.80 


0.006 







010 


0.80 


- 1.00 


0.006 







010 


1.00 


- 2.00 


0.006 







010 


2.00 


- 3.00 


0.005 







010 


3.00 


- 4.00 


0.003 







005 


4.00 


- 5.00 


0.003 







005 


5.00 


- 6.00 


0.003 







005 


6.00 


- 8.00 


0.0017 







003 


8.00 


- 10.00 


0.0017 







0017 


10.00 


- 12.00 


0.0017 







0017 



the CMDs of the association. First, the CMDs of both 
association and field are divided in boxes of the size of 
AF435H^ = 0.2 and A(F435T^ - F8UW) = 0.3. The in- 
completeness correction is taken into account by dividing 
the field and association CMDs in magnitude-color bins 
and then adding on each bin having Nth stars, (1 — A) x Nth 
objects, where A is the smallest of the FA35W and F8MW 
completeness factors. At each CCD area in both the as- 
sociation and in the field the appropriate completeness 
correction as calculated from artificial star experiments is 
applied. Then the closest cluster star is subtracted in each 
box of the Nil CMD for every field star. This procedure 
was already used in a number of papers where more de- 



tail c an be found (jChiosi et al 
2003). 



2006aHChiosi fc Vallenari 



4.3. Determination of the interstellar extinction in 
N 11 

In order to study the stellar population in the association, 
it is necessary to derive a detailed map of the interstellar 
extinction which is expected to be highly variable. First, 
we select the main sequence stars brighter than F435W 
~ 18. The comparison between the CMDs of the field 
population and those of Nil shows that these stars have 
a high probability to belong to the association, since they 
are brighter than the youngest turnoff of the field popu- 
lation (see Figs. [5] and [3]) . About 1000 stars are found in 
this magnitude range. Then, assuming that the spread in 
color of the stars at a given magnitude is due to inter- 
stellar extinction only, their color (F435W^ — F814W^) is 
compared with the intrinsic color expected for main se- 
quence stars having an age of 5 Myr (see discussion in 
Sect. 6), which allows the determination of the color ex- 
cess E(i^435W^ - i^814W^) for each star. We assume the 
extinction law Ap^^^w /^fsuw by Sirianni et al (2005). 
Then the E(i^435M^ - FSUW) values are interpolated us- 
ing the package GRIDDATA in IDL as a function of the 



coordinates. Finally, the extinction of each star in the field 
is calculated from the interpolated map. 

4.4. Statistical analysis of the stellar spatial 
distribution 

We present the methods used to study the degree of clus- 
tering of the stars and derive the spatial scale of the forma- 
tion. First, we define the two-point correlation function, 
second we describe the minimal spanning tree method. 

4.4.1. Two-Point correlation function 

The probability 1 -|- ^(r) of finding a neighbor in a shell 
element at a distance r from an object of the sample is 



N 



l+ar) = l/iNn)J2Mr). 



(1) 



i=l 



where ^(r) is the two-point correlation function, ni(r) 
is the number density of objects found in an annulus cen- 
tered on the i-th object and having a radius between r and 
r-f-dr, N is the total number of ob jects, and finally, n is 
the average number (|Peeblesl . [l980l) . A Monte Carlo algo- 
rithm is used to derive the area covered by the data when 
the annulus extends outside the studied region. Using the 
above definitions, a random distribution of stars will pro- 
duce a fiat correlation, with ^(r) 0. A peaked ^(r) at 
small radii indicates a positive correlation, and the full 
width half maximum of the peak itself represents the spa- 
tial scale of the clustering. The distribution of the near- 
est neighbor distance might reflect the structure of the 
interstellar medium from which the stars formed. The ab- 
solute value of ^(r) is a measure of the concentration of 
the objects at a given dist ance r relative to the average 
distribution (jPeebled . [l980[) . 



4.4.2. Minimal Spanning Tree 



The Minimal Spanning Tree (MST) is a powerful 
method used in literature for a statistical analysis 
of the spatial distribution of the stars in a cluster 



(see 



2006 



Cartwright fc Whitworthl . 120041; ISchmeia fc Klessenl 



for a detailed discussion). The MST is the unique 
network of straight lines joining a set of points (vertexes), 
so that it minimizes the total length of all the lines (called 
edges) without creating clo sed loops. In this pap er, we 
make use of the routi ne ms^ ( Havmond et al.lll984l ). based 
on Prim's algorithm ( Priml . 1957^ to construct the MST. 
Starting from an arbitrarily chosen first point, an edge is 
created joining that point to its nearest neighbor. The tree 
is then calculated constructing the shortest link between 
one of its nodes and the unconnected points, until all the 
points are connected. Once a MST is calculated, its mean 
edge length m is directly defined. The m is dependent on 
the total number N of stars in the cluster, in the sense 
that as N increases, sho rter edges are created and to de- 



creases. For this reason, ICartwright fc Whitworthl (|2004f ) 
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normalize m to {NAY^'^ /(N— 1) when comparing clusters 
that have a different obj e ct nu mber N and/or areas A. 
Cartwright fc WhitworthI ( 2004 ) and ISchmeia fc Klessen 



(l200d) demonstrate that the structure of a cluster can be 
derived comparing the m of the MST describing it, to the 
normalized correlation length s defined as the mean sep- 
aration between the stars normalized to the radius of the 
region. 

This defines the parameter Q as Q = ™ . 



Cartwright fc WhitworthI ()2004D show that s decreases 



more quickly than fh as the degree of the central concen- 
tration becomes more severe (indicating the presence of a 
radial gradient), while m decreases more quickly than s as 
the degree of subclustering is getting higher. Q can be used 
to distinguish between subclustering ( Q < 0.8) and the 
presence of a large scale radial density gradient (Q > 0.8). 
It should be remembered that Q cannot give any infor- 
mation about the fractal nature of the sub-clustering, i.e. 
it cannot say whether the subclustering is hierarchically 
self-similar. 




5. Results: the SFR history of the field population Fig. 4. Extinction map in E(F435W-F814W). 



Fig. [3] presents the CMD of the field population where 
isochrones of different ages are superposed. The youngest 
field population shows a turnoff at about F435W ~ 17—19, 
and the stars burning He in the core are located at F435W 
~ 20. The CMD suggests the presence of an intermediate 
age and young population. 

The detailed determination of the star formation rate 
is presented in Fig. [3] (lower panel) and reveals the pres- 
ence of a small component of a very old population 
(10—12 Gyr), an intermediate age component in the range 
2—6 Gyr, and several bursts of a younger population, up 
to 0.4 Gyr. Th is result agree s with pre vious dctcrmina- 
tions bv [Dolphin et all (|200ll) . lHarr"is&: Zaritskvl (|2004l ): 



Javiel et al.l (|2005l) . 



6. Results: the interstellar extinction 

The map of the interstellar extinction in E(F435VF — 
F814W^) is presented in Fig. |4l The interstellar extinc- 
tion E(F435W-F814W) is found to be in the range of 
0.2—0.5. This corre sponds with Afa^sw ~ 0.95-0.38 us- 
ing the relations bv lSirianni et al. ( 2005 ) and agrees with 
the value s derived by ni e ans o f the spectroscopy of bright 
stars by iMokiem et al. ( 2007t ). The main source of un- 
certainty on our determination of the extinction is due 
to the presence of an age spread inside the young main 
sequence stars. We derive the extinction using main se- 
quence isochrone with an age of 5 Myr, while ages as young 
as 2 Myr are found for the blue population (see discussion 
in the following Section). At F435W=12.5, the expected 
MS color difference between an age of 2 Myr and 5 Myr 
is of 0.07 mag, while at F435W=14 is 0.01 mag. These 
values can be regarded as an estimate of the uncertainty 
on'E{FA35W - F8UW). 




F435W-F814W 



Fig. 5. CMD of the N 11 region covered by the HST/ACS 
photometry. The lines drawn in the CMD mark the regions 
redder than the MS where PMSs are clearly separated 
from the field population. These lines indicate only a part 
of the region in the diagram where PMSs are located. 



7. Results: Ages and stellar content of the 
associations inside N 11 

Here we analyze the CMDs of the most relevant associa- 
tions inside N 11 with the aim of deriving the age of the 
stellar population. We assume a mean value of Z=0.19 for 
the metal content of the association, in agreement with 
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NGC1760 
R^l.larcmin 




Age- 2, 5. 10 Myr 
A„^^w= 0,7 




(F435W-Fai4W) 



{F435W-F814W) 



HD32228 
R-0.55arcmin 



Age= 3, 5, 10 Myr 
A,,,„.^= 0.5 




(F435W-F814W) 



(F435W-F814W) 



Fig. 6. CMDs of NGC 1760, HD 32228 in LH 9 (left panels). CMDs of the same objects when the field population is 
subtracted (right panels). The bars show the photometric errors on the color. PMS isochrones by Siess et al. (2000) 
are superimposed on the data. The arrow indicates the extinction vector. 



Pagel fc Tautvaisiend (|l999l ). IPirsch et all (|200d ). 



Figure [5] present the CMD of the Nil region covered 
by the HST/ACS photometry. A well-defined upper main 
sequence of young blue stars is present. Below the turnoff, 
moving to fainter magnitudes, the main sequence in the 
CMDs of the associations and of the field population is in- 
creasingly densely populated. A striking difference at faint 
magnitudes between field and association CMDs is evident 
in the region redder than the MS, where a well-populated 
sequence almost parallel to the MS itself is found in the 
associations, but is completely absent in the field region 
(see Figs Island El . These stars fit the location of the PMSs 
previously discovered in other MC associations very well 
i|Gouliermis et all . I2007ll2006l ). 



Suspicion might arise that the stars we identify as PMSs 
are instead highly reddened MS objects. However, the fact 
that the interstellar extinction vector in the passbands 
F814W-F435W is almost parallel to the PMS sequence(see 
Figs. [6l [71 [To]) means that the interstellar extinction has 
only a small effect on the age of the PMSs. 

7.1. The stellar population in LH 9 

The most relevant associations/clusters in LH 9 in the 
region observed by HST are NGC 1760, NGC 1761, 
HD 32228, BCDSP 1, BSDL 270. 

Figures [HI and [3 present the CMDs of these objects 
before and after the field star subtraction. The PMS pop- 
ulation is clearly visible in all the fields. The CMDs after 
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BCDSP1 
R=0.65arcmin 



AgR= 2. 5, 10 Myr 




(F435W-FB1^W) 



(F435W-F814W) 



BSDL270 
R=0.46arcmin 



Age= 2, 5. 10 Myr 




(1'435W - F814W) 



fF435W-F814ff) 



Fig. 7. CMDs of BCDSP 1, BSDL 270 in LH 9 (left panels). CMDs of the same objects when the field population is 
subtracted (right panels). The bars show the photometric errors on the color. PMS isochrones by Siess et al. (2000) 
are superimposed on the data. The arrow indicates the extinction vector. 



the field star subtracti on are compared w ith pre-main se- 
quence isochrones by ISiess et ah 



2000 ^ converted frc 



the Johnson-Cousins pass-bands t o the HST/ACS pass - 
bands using the transformations bv lSirianni et al" ( 2005 ). 
Below we discuss the age determination of each object. As 
a general remark, we remind that it is difficult to derive 
a precise age determination of the PMSs on the basis of 
the CMDs only. Several effects such as photometric er- 
rors, binarity, differential extinction across the field and 
age spread in the populat ion cannot be disentangl ed (see 



for a detailed discussion Hil lenbrand et al. , 12008'). This 



holds for all the objects studied here. This problem will 
be discussed in the following sections. 



NGC 1760. Comparing the location of the PMSs with 
isochrones, we derive a tentative estimate of the age in the 
range of 2—10 Myr. (see Fig.[6|). More stringent constraints 
can be derived on the basis of the luminosity function. The 
observational LP after incompleteness correction, field star 
subtraction and extinction correction is compared with 



the theoretical one expected from the ISiess et al. (|200fl[ ) 



PMS isochrones. We make use of an authomatic proce- 
dure to find the best fitting solution. The LP of the PMSs 
(see Fig. [H]) is suggesting a prolonged formation from 2 to 
10 Myr (for single stars) with a IMF slope P = -3.0 ± 0.3 
(when the Salpeter value is P = —2.35). In the best-fit 
solution, the majority of the stars (about 75% of the ob- 
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Fig. 8. Observational LF of NGC 1760 corrected for in- 
terstellar extinction, completeness, and field star subtrac- 
tion (solid line) is compared with theoretical LF for a 
prolonged star formation in the age range 2—10 Myr as 
derived using Siess et al. (2000) isochrones (dashed line) 
for stars fainter than F435W ^ 17. 



jects brighter than F435W ~ 24.5) are however older than 
6 Myr, while the youngest population (2 Myr) can account 
for less than 1% of the population. The mass range of the 
PMSs is from 1.5 to 3.5 Mq and 1.38 to 1.7 Mq for 2 and 
10 Myr, respectively. The lower mass limit corresponds to 
F435W ^ 24 . 0. Ou r value of the IMF slope agrees with 
Parker et all (|l992[ ) who find T = -2.6 ± 0.1 for masses 
higher than 8-9 M©. 

The question at hand is how reliable the age and age 
spreads are which we find from the CMD and LF analy- 
sis. With regard to the age spread, photometric errors and 
binarity can complicate the analysis. Binarity can be dis- 
tinguished from age spread only when empirical lu minos i- 
ties are known with high accuracy (Hillenbrandl . |2009| ). 
For single stars, taking into account photometric errors, 
we can say that the percentage of stars older than 6 Myrs 
is ranging from 65% to 83%. On the basis of this dis- 
cussion, it would be unrealistic to reach a more precise 
conclusion concerning the age distribution. Multiplicity 
surveys found a relatively high binary frequency ( >60%) 
for masses higher tha n solar in young associations in the 
Galaxy (among others iDuquennov fc Mavor , 199lh . If we 
impose that 60% of stars are in binary systems, we do not 
derive a significantly different age spread. 

In addition to these effects, it is well known that sys- 
tematic trends are found in the literature using various 
sets of pre-main sequence tracks. This holds in particular 
at sub-solar masses, while at solar masses the agreement 
is better. In addition, an age-with-mass trend is found, 
i.e. higher stellar masses are predicted to be older than 
lower mass stars in the same objects using the same set 
of stellar tracks. All this can easily account for 20% to 
100% of the age difference when comparing the ages we 
derive with values taken from various literature sources 

2008t iHillenbrand 



Fig. 9. Observational LF of HD 32228 (solid hue) is com- 
pared with theoretical LF for single stars derived using 
Siess et al. (2000) isochrones (dashed line) for stars fainter 
than F435W 17. (see text for details) 



a remarkably good agreement of our PMS ages with the 
literature spectroscopic determination of t he ages of blue 
stars. OB stars are found in NGC 1760 bv iMokiem et al ' 



(|2007[ ). This sets a limit to the age at 7.0 ±1 Myr. These 
authors find that the age range is broad in the whole as- 
sociation, going from about 7 to 1—2 Myr. 

HD 32228. A high concentration of PMSs is found in 
HD 32228 located at the periphery of the association LH 9. 
Comparing the location of the PMSs with isochrones we 
find that the age of the PMS is in the range of 3—10 
Myr (see Fig. [6]). If we assume that no binary stars are 
present, the LF of HD 32228 PMS after field star sub- 
traction and extinction correction suggests that the dom- 
inant population has an age of 5 ±1 Myr and a IMF slope 
r = -2.0 ± 0.2 (see Fig. ^. The mass range of the PMSs 
reaches from 1.5 to 2.5 Mq when F435W < 24.0. However, 
a small age spread cannot be ruled out: solutions includ- 
ing, in addition to the main population of 5 Myr, about 
25% of the stars are as old as 7 Myr, and a small fraction 
younger than 5 Myr cannot be ruled out. Assuming that 
about 60% of the stars are in binaries, we derive a slightly 
older age of the dominant population (6 ± 1 Myr), with 
a Salpeter IMF. The age of the cluster HD 32228 derived 
from spectroscopic of bri ght stars is of about 3—4 Myr 
(jWalborn fc Parkeii Il992f ) based on the presence of a WC 
star. 

BSDL 270 and BCDSP 1. A small population of 
PMSs is found in these two small objects located at the 
periphery of LH 9. Comparing the location of the PMSs 
with isochrones in BSDL 270 and BCDSP 1, we find that a 
reasonable estimate of the age is in the range of 5—10 Myr 
(see Fig. [7]). Due to the scarcity of the PMS number, the 
LF cannot reliably give information about the age of the 
stars. Several B stars are found in these clusters, and one 



06.5V (Nll-065) is in BCDSPl (jMokiem et all 120071 ). 



(see among others IHillenbrand et al 



200i). Taking these uncertainties into account, we find 



Nll-065 has a location on the HRD compatible wi th an 
age of about 1 Myr. However, iMokiem et al.l (|2007l) find 
that this object has high He abundance. Using the ap- 



10 



A. Vallenari et al.: Pre-Main Sequence stars in N 11 in the LMC 



NGCir63 
R=1 .59arcmin 



AgE= 2. 5, 10 Myr 
Af,^™= 0,8 




1 S 

(F435W-FB14W) 



(F435W-F814W) 



Fig. 10. CMD of NGC 1763 (LH 10) (left panel). CMDs of the same objects when the field population is subtracted 
(right panel). The bars show the photometric errors on the color. PMS isochrones by Siess et al. (2000) are superimposed 
on the data. The arrow indicates the extinction vector. 



propriate isochrones for the surface He abundance, they 
derive an age of 6 Myr. It is quite remarkable that in the 
whole region inside the uncertainties there is a good agree- 
ment between the ages of young blue stars and of the red 
PMSs. 



7.2. The stellar content of LH 10 

The CMD of NGC 1763 in LH 10 before and after the field 
population is subtracted is shown in Fig. [TOl The CMDs 
refers to a small region near the center of a radius R ^ 1.6'. 
PMSs are clearly identified in the CMD. The LF for single 
stars is presented in Fig. [TT] after incompleteness correc- 
tion, field star subtraction and extinction correction. The 
comparison of the observ a.tional mod e ls wit h the theoret- 
ical one calculated using ISiess et al. models sug- 
gests that the dominant population has an age of 2—3 
Myr, using a IMF slope F = -2.20 ± 0.2. The PMS mass 
range is 1.5—3.5 M0. As we discussed in the previous sec- 
tion, it is difficult to completely discard solutions having a 
small percentage of older objects. This result agrees with 
Parker et al.l (jl992l ) who find that the slope of the initial 
mass function of the young blue MS stars in LH 10 is 
F = -2.1 ± 0.1, significantly flatter than that of LH 9. If 
we account for 60% of binaries, we derive an age of 3 ±1 
Myr, wfth a steeper I MF (F = -2.80ib0.2). Us i ng br ight 
blue star spectroscopy, iHevdari-Malaveri et al. I(l2000h es 



timate that LH 10 in NllB is 3±1 Myr old. lMokiem et al 



(|2007l ) derive a broad age spread in LH 10 from 1 to 6 Myr. 
But the majority of the stars are younger than 4.5 Myr, 
and only one object (Nl 1-087) in the cluster core is older 
than this age. These authors conclude that the mean age 
is 3±1 Myr. Herbig Ae/Be stars are found inside the as- 



sociation, together with a methanol maser (jHatano et al.l . 
.2006.) . 

7.3. The stellar population in LH 13 

Our HST data in LH 13 include NGC 1769 and BSDL 324. 
The CMDs of the association before and after the field 
population is removed are shown in Fig. 1121 The compar- 
ison of the location of PMS on the CMD of NGC 1769 
with isochrones suggests an age in the range of 2— 5 Myr. 
Analyzing the PMS LF (see Fig. [13]) we find that the dom- 
inant population (about 85% of the stars) have an age of 
2—3 Myr with a Salpeter IMF when only single stars are 




Fig. 11. Observational LF of NGC 1763 (solid line) is 
compared with theoretical LF (single stars) derived using 
Siess et al.(2000) isochrones (dashed line) for stars fainter 
than F435W-17. 
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Fig. 12. CMDs of NGC 1769 and BSDL 324 (in LH 13) (left panels). CMDs of the same objects when the field 
population is subtracted (right panels). The bars show the photometric errors on the color. PMS isochrones by Siess 
et al. (2000) are superimposed on the data. 



considered. The PMS mass range is 1.5—3.5 M0. Including 
60% of binaries, we derive an age of 3—4 Myr, with a steep 
IMF slope F — — 2.8±0.1. No information about the age of 
the PMS in BSDL 324 can be derived from the CMD, due 
to the small number of stars detected. Spectroscopic deter- 
mination of the age of OB stars in LH 13 are in the range of 
3—5 Myr (jHevdari-Malaveri et al. ■ l2000h . A small concen- 



tration of Herbig Ae/Be stars is found here (|Hatano et al 
2n06h . 



7.4. YSOs from near-IR archive data: data selection 
and object distribution 

In this section the detection of YSOs on the basis of 
IR Spitzer archive data is presented. The comparison of 
the observational color-col or plots in the i nfrare d pass- 
bands with the models by iRobitaille et al. 1 dlooi) allows 
us to select candidate YSOs of Stage O-I (presenting sig- 
nificant infalling envelope , and having ages of < 0.1 
Myr) and Stage II (showing a thick disk source, hav- 
ing 0.1 <age< a few Myr) and finally. Stage III YSOs 
sources (with an optically thin disk or no disk ) . Her e 
we adopt the classification by IRobitaille et al. I (|2006l) . 
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which is analogous to the traditional Class scheme (jLada , 
19871) . On the basis o f 200,000 YSQ models computed 
at ten viewing angles, ^ Robitaille et al. 1 (l2006l) point out 
that the IRAC color-color plot ([5.8] [8.0])-([3.6] [4.5]) con- 
tains a large region where the fraction of Stage 0-I/all 
YSOs is close to 100%, a small intermediate region where 
the fraction of Stage 11/ all sources is higher than 80%. 
IRAC-MIPS ([8.0][24])-([3.6][5.8]) are effective in sepa- 
rating stars with no circumstellar material from Stage 
O-I and Stage II objects. Here we restrict ourselves to 
study Stage O-I and Stage II objects. In Figs[T4] and [15] 
the SAGE archive data in the N 11 region are plotted 
in the IRAC color-color plot ([5.8] [8.0])-([3.6] [4.5]) and 
IRAC-MIPS color-color ([8.0][24])-([3.6][5.8], respectively 
to identify the YSOs. Finally YSOs candidates selected 
on the basis of previous plots are examined in the [8.0]- 
([3. 6]- [8.0]) color-m a gnitu de diagram and compared with 
the IWhitnev et"al] \2004 i models (see Fig. [TH). About 
30 objects are identified as YSO candidates. Background 
galaxies and foreground Milky Way stars are known to 
contaminate the IRAC and MIPS CMDs and color-color 
plot. In order to derive the expected contribution in the 
studied region, we select a control field located outside 
the star forming region, namely at a(J2000)=4h 58m 40s, 
5(J2000) -66° 34/ 47", which has an area comparable to 
the area of the whole Nil field. Using the same diagnostic 
plots we find six objects in the control field redder than 
([3.6]-[8.0])'- 1 in the [8.0]-([3.6]-[8.0]) color-magnitude di- 
agram and only 3 (over 30 in total) in the region oc- 
cupied by YSO candidates in the ([5.8][8.0])-([3.6][4.5]) 
color-color plot. Only 10% of the YSO candidates having 
[8.0] < 14 are expected to be background objects. The per- 
ce ntage of background g alaxies agrees with the statistics 
bv lMeixner et all (|2006l ). 

The YSO spatial distribution is shown in Fig. (TH] type 
I and II objects having ages from 0.1 to 1 Myr are found 
at the same location as the PMSs, indicating that sev- 
eral generations of stars are present in the region. The 
youngest candidates (type I) are found in NGC 1760 
(LH 9), north of BCDSPl (LH 9) and NGC 1769 (LH 13). 



7.5. Toward a coherent picture: the PMS, YSOs, OBs, 
Herbig Ae/Be and their age distribution 

In Fig. [T7] we compare the spatial location of the PMSs 
with the location of Bica et al clusters and associa- 
tions in the area and with the distribution of t he OB 

20061 ). 



stars and Herbig Ae/Be candidates ([Hatano et al. 
Summarising the results of the previous sections, the ages 
derived in literature from spectroscopy of blue stars and 
from PMS analysis are quite consistent in the whole area. 
The analysis of HST data imphes that PMSs in LH 9 
South the majority of the stars are slightly older than the 
remaining part of N 11, while LH 9 North, LH lO(center) 
and LH 13 seem to be younger (with ages of 5, 2, 2 Myr). A 
large age spread from 1 to 6 Myr is found from the spec- 
troscopy of blue stars in these regions and is confirmed 




Fig. 13. Observational LF of NGC 1769 (solid line) is 
compared with theoretical LF derived using Siess et al. 
(2000) isochrones (dashed line) for stars fainter than 
F435W - 17. 
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- - Robitaille et al 2006 Stage II 




- - Robitaille et al 2006 Stage 1 
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Fig. 14. [5.8][8.0])-([3.6][4.5] color-color plot of the N 11 
field. Red identifies YSO candidates. Error bars show the 
nominal errors on the photometry as given by the SAGE 
catalog. The blue line outlines the region wh ere Stage O-I 
objects are found on the basis of iRobitaille et al.l ( 20061 ) 
photometric models. In the region redder than [5.8] [8.0] 
~0 between the green and the blue line Stage Il/Stage HI 
objects can be detected. 

by our analysis. The presence of YSOs in LH 9 South, 
LH 13 (type I) and LH 9 (North), LH 10 (Type II) sug- 
gests that the star formation in the region is still going on. 
The most conspicuous concentrations of PMSs correspond 
to the location of OB stars. Herbig Ae/Be stars are found 
in all these locations, confirming the presence of a popu- 
lation of ages in the range of 1—3 Myr. All the indicators 
prese nt a coherent picture of the star formation in the re- 
gion. iMokje^^eT^l] (2007) estimate that the time needed 
for a supernova shock to cross the distance between LH 9 
and LH 10 is of the order of 0.1 Myr. Since the lifetime 
of a massive star is of the order of 3 Myr, the age dif- 
ference agrees with the idea that LH 9 has triggered the 
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Fig. 15. IRAC-MIPS color-color plot in the pass-bands 
([8.0][24])-([3.6][5.8]). The short dashed red line separates 
th e region where S t age O- I objects are found, on the basis 
of iRobitaille et al I (l2006l) photometric models. The long- 
dashed green line shows the region where Stage II objects 
can be detected. 



• Young Stellar Object Candidates 
Wliltney et al 2004 YSO 0-1 




56m48s 
a(J2000) 



Fig. 17. Bica clusters and associations are overplotted on 
the distribution of PMS (dark magenta dots) , the concen- 
tration of OB stars (light blue circles) and Herbig Ae/Be 
stars (orange triangles) (from Hatano et al 2006). Red el- 
lipses represent Bica et al. object dimensions (see Fig. 1 
for details). The blue contours show the observed area. 




Fig. 16. [8 . 0] - ( [ 3 . 6] - [8 . 0] ) color-mag nitude diagram com- 
pared with the I Whitney et all (2004) models. The objects 
redder than the line are likely to be YSOs of the Stage 
O-I. 



formation of a first star generation in LH 10. This cannot 
be excluded on the basis of our discussion, even if more 
generations of stars are found in the region. 

We compare the location of PMSs and YSOs with 
emission (Mac Low et al 1998) and CO emissions (Israel 
et al 2003) (see Fig. [HI). We find that PMSs and YSOs 
are located near the maxima of the emission and at 
the border of the CO clouds in NllB and NllC. PMSs 
are found at the border of the cavity surrounding LH 9. 



8. Results: clustering of PMS stars 

In this section we discuss the clustering degree of the 
PMSs, first using the two-point correlation function and 
then the MST method. Only the stars selected as PMSs 
on the basis of the HST/ACS photometry are considered. 

The correlation function for the PMSs shows a well 
defined peak, while in the case of the whole star sam- 
ple it is considerably flatter, implying a weaker corre- 
lation (see Fig. [T9)) . The maximum value of the cor- 
relation function of the PMSs is comparable to the 
peak values found in Galacti c star formi n g reg ions such 
as Rosetta Molecular cl oud (|Li fc Smithl . boosh or M16 



(jlndebetouw et al.l . 120071) and reflects the strong degree of 



clustering of the population. 

The full width half maximum of the peak for PMSs is 
of ^ 10 pc, of the order of the typical size of molecular 
clouds in the SMC and LMC w hich is going from about 



10 to 100 pc (|lsrael et al.l . l2003h . This reflect the structure 



of the interstellar medium from which the stars formed. 
This result may however be affected by a potential bias 
due to the fact that the ACS pointings are not covering 
the whole region, but are mainly centered on well-known 
star forming regions, where higher concentrations of young 
stars are expected. 

In Fig. [20] we present the MST graphs of some of the 
clusters and associations in the region. In Table [3| we give 
the values of the Q parameter, calculated selecting the 
stars inside the elliptical areas described by the Bica et 
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al. cluster parameters. Although in princ iple Q is not 



expec ted to depend on the selected area, iBastian et al, 
( 2009f ) found a degeneracy in the sense that highly elon- 
gated structure like these show a lower Q value. The val- 
ues hiJT^\e_^\jnclude the ellipticity correction derived 
by iBastian et al. ( 20091 ). To avoid the uncertainties due 
to the shape, we calculate Q using the quadrilateral ar- 
eas which correspond to the observational CCD fields. The 
results are presented in Table S] for all the stars. No signif- 
icant changes are found due to the area selection. To ver- 
ify whether blue bright stars and fainter red PMSs have 
a similar clustering behavior, we calculate Q for PMSs 
fainter and brighter than F435W ~ 20. Tables M and [6] 
respectively show the results. We use quadrilateral areas. 
Finally we derive Q for the OB star candidates taken from 
the Hatano et al. catalog in an area of 0.24 sq. deg. in the 
centre of the Nil star forming region. The corresponding 
graph is presented in Fig. [20l The Q parameter is derived 
for the whole area, without subdivision in smaller fields, 
due to the small number of OB candidates. Figure [H] 
shows the trend of Q with the magnitude of the popu- 
lation. Q is calculated for all the stars brighter than a 
given magnitude. The main results can be summarized as 
follows: 

— In all the fields related to Nil, the values of Q 
< 0.8 are typical of clustered substructures. Smaller 
objects have slightly scattered Q values in the 
range of 0.62— 0.74. For compar ison, we remind that 
ICartwright fc WhitworthI (|2004h find that a value of 
Q ~ 0.56, 0.66, 0.70 corresponds to a projected fractal 
model with dimensions of 1.8,2.2, 2.4 respectively. 

Table 3. Values of Q for the all the stars inside the cluster 
areas given by the Bica et al parameters. 



Table 4. Values of Q for all the stars in quadrilateral areas 
corresponding to the CCD fields (see text for details). 



Object 


m 


s 


Q 


NGC 1769 


0.86 


1.33 


0.65 


NGC 1763 


0.73 


0.98 


0.75 


NGC 1761 


0.77 


1.18 


0.66 


NGC 1760 


0.73 


1.20 


0.56 


HD 32228 


0.73 


0.93 


0.79 


BSDL 324 


0.80 


1.11 


0.73 


BSDL 270 


2.45 


3.64 


0.67 


BCDSP 1 


0.78 


1.19 


0.66 



OB candidate distribution has Q= 0.85. This indi- 
cates a diffuse population. Two main concentrations 
are found however, associated to HD 32228 (and in 
general to NGC 1761) and NGC 1763. In these ob- 
jects, the positions of the concentration of OB stars 
are very close to the main concentration of PMS stars. 
Inside Nil no significant difference is found in the 
clustering degree of brighter blue main sequence stars 
(F435W^ < 20) and fainter red stars (F435W > 20) 
which suggests that they are formed in the same pro- 
cess. A closer look reveals that the regions present 



Object 


m 


s 


Q 


NGC 1763 


0.68 


0.98 


0.70 


NGC 1769 


0.70 


1.13 


0.62 


LH 9 East 


0.71 


0.97 


0.74 


(NGC 1761 E) 








LH 9 North 


0.71 


0.97 


0.74 


(NGC 1761 W,HD 32228, 








BCDSP 1,BSDL 270) 








LH 9 South 


0.71 


0.97 


0.74 


(NGC 1760) 








Field stars 


0.76 


1.12 


0.67 



Table 5. Values of Q for PMS (fainter than F435W - 20 
) in quadrilateral areas corresponding to the CCD regions. 



Object 


m 


s 


Q 


NGC 1763 


0.70 


1.00 


0.70 


NGC 1769 


0.71 


1.15 


0.62 


LH 9 East 


0.77 


1.11 


0.70 


(NGC 1761 E) 








LH 9 North 


0.75 


1.04 


0.72 


(NGC 1761 W,HD 32228, 








BCDSP 1,BSDL 270) 








LH 9 South 


0.75 


1.12 


0.67 


(NGC 1760) 









a different behavior concerning the clustering degree 
as a function of magnitude (see Fig. [2T|) . While in 
NGC 1763, LH 9 East and North bright stars are more 
concentrated, in NGC 1769 ad LH 9 South are quite 
dispersed. PMS fainter than F435W ^ 20 do not show 
any trend with regard to the magnitude. 
Since young clusters are expected to evolve from hi- 
erarchical configurati ons to more concentrated ones 
( Schmeia et al. . 20091 ) . the fact that there is a range 
of structures with a scatter in Q may indicate that the 
objects are at different evolutionary stages. 
Q is found to have no dependence on the age of the 
clusters, even if the oldest objects HD 32228 and in 



Table 6. Values of Q for BMS (F435W<20) stars (quadri- 
lateral areas). 



Object 


m 


s 


Q 


NGC 1763 


0.63 


0.93 


0.68 


NGC 1769 


0.67 


1.12 


0.58 


LH 9 East 


0.76 


1.09 


0.69 


(NGC 1761 E) 








LH 9 North 


0.70 


0.97 


0.72 


(NGC 1761 W,HD 32228, 








BCDSP 1,BSDL 270) 








LH 9 South 


0.76 


1.16 


0.65 


(NGC 1760) 
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Fig. 18. Ha emission map (inverted intensity scale) by McLow et al (1998). Large boxes show the ACS/WFC fields. 
Green contours indicate the CO clouds (Israel et al 2003), small red dots show the PMS candidates selected on the 
basis of ACS photometry, large green squares show the location of YSOs type II selected using IR Spitzer photometry, 
large blue dots repre sent the location of Y SOs type I. Finally red circles (marked by HCl to IIC6) show the Herbig 
Ae/Be candidates bv lHatano et al. ( 20061 ). 




50 100 
radius (pc) 



general LH 9 North have a slightly higher Q. It should 
be remembered that this comparison is complicated by 
the presence of a large age dispersion detected in the 
fields. 

This is at odds with what we find in the star form- 
ing region in our Galaxy, where Q correlates with the 
age. ISchmeia et al. I (l2009l ) find similar behavior in the 
region of NGC 346 in the SMC and suggest that this 
can be explained as the footprints of inhomogeneity in 
the original turbulent cloud. In fact Q is expected to 
depend on a variety of additional factors like the initial 
conditions or turbulent energy in the clouds. 
The value of Q itself in the field population is rather 
low, around 0.67, meaning that the field still has sub- 
structures due to the original turbulent environment 
or showing the result of random motion configurations. 
The value we f i nd is compatible with Q=0.7 derived by 
Bastian et al. ( 2009f ) as a reference value for the field 



Fig. 19. Two-point correlation function of the PMSs 
(solid line) and of the whole sample of stars (dashed line) . 



of the LMC. 
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Fig. 20. The MST of the regions of NGC 1769,NGC 1763,NGC 1761, NGC 1760, HD 32228, BSDL 324, BSDL 270 
and BCDSP 1 and for OB star candidates. 



9. Conclusions 

Young stars are well known to trigger star formation on 
a local scale by injecting energy in the surrounding inter- 
stellar medium. Nil region in the LMC is often presented 
as one of the most promising candidates of a triggered 
star formation. In this paper we discussed the star for- 
mation process in N 11. Wc reported on the discovery of 
PMS and YSO candidates associated with Nil from HST 
ACS/WFC photometry and archive Spitzer data. While 



HST observations can give information about faint, ex- 
posed prc-main sequence candidates, IR data allowed us 
to detect embedded young stellar objects. The main con- 
clusions are as follows: 

— A large population of PMSs is found in this region. 
The comparison with the isochrones shows that PMSs 
have masses from 1.3 Mq to 2.0 M0, depending on 
their ages. YSOs (type I and II) having ages < 1 Myr 
are found at the same location as the PMSs. 
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— The analysis of HST data impHes that PMSs in LH 9 
are slightly older than in the other regions of N 11. 
A prolonged star formation is detected in LH 9 South 
from 2 to 10 Myr. However, the majority of the stars 
are older than 6 Myr, while LH 9 North, LH 10(cen- 
ter) and LH 13 seem to be younger (having ages of 
5, 2, 2 Myr respectively). A large age spread from 1 
to 6 Myr is found from the spectroscopy of blue stars 
in these regions. Remarkably, the present age deter- 
minations of the PMSs agree (inside the errors) with 
the ages derived from literature spectroscopy of blue 
bright stars. 

— The presence of YSOs in LH 9 South, LH 13 (type I) 
and LH 9 (North), LH 10 (Type II) suggests that the 
star formation in the region is still going on. The most 
conspicuous concentrations of PMSs correspond to the 
location of OB stars. Concerning Herbig Ae/Be stars, 
they are found in all these locations, confirming the 
presence of a population with ages in the range of 1—3 
Myr. All the indicators present a coherent picture of 
the star formation in the region. 

— PMSs present a higher degree of clustering in compari- 
son with the whole star sample. The distribution of the 
nearest neighbor distance reflects the structure of the 
interstellar medium. The full width half maximum of 
the peak of the correlation function for PMSs is ^ 10 
pc, which is the order of the typical size of molecu- 
lar clouds in the SMC and LMC, which extends from 
about 10 to 100 pc. 

— In all the fields related to Nil, the values of Q < 0.8 
are typical of the products of hierarchical star forma- 
tion. Clusters are expected to evolve from hierarchical 
structures to more concentrated ones. Different clus- 
ters/associations show slightly different values of Q, 
suggesting that they might be in slightly different evo- 
lutionary stages. Inside Nil, no significant difference 
is found in the clustering degree of brighter blue main 
sequence stars (F435VF < 20) and fainter red stars 
(F435W > 20), suggesting that they might be part of 
the same formation process. 

— Q is found to have no dependence on the age of the 
clusters, even if the oldest objects HD 32228 and in 
general LH 9 North have a slightly higher Q. These 
results disagree with the behavior found in our Galaxy. 
These results are quite difficult to interpret. They can 
be explained as the footprints of inhomogeneity in the 
original turbulent cloud. It should be kept in mind 
that this comparison is complicated by the presence of 
a large age dispersion detected in the fields. 
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